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ABSTRACT 

We measure reddenings to H II regions in NGC 2903, NGC 1512, M51, NGC 4449 
and NGC 6946 from Hubble Space Telescope Paa and Ha images. Extinctions range 
from Ay ~ 5 ~ depending upon the galaxy. For the galaxies with HST images in 
both lines, NGC 2903, NGC 1512 and M51, the Paa and Ha emission are almost 
identical in morphology which implies that little emission from bright H II regions 
is hidden from view by regions of comparatively high extinction. The scatter in the 
measured extinctions in each galaxy is only ±0.5 mag. 

We compare the reddenings we measure in five galaxies using the Paa to Ha ratios 
to those measured previously from the Balmer decrement in the LMC and as a function 
of radius in MlOl and M51. We find that luminosity weighted mean extinctions of 
these ensembles of H II regions are correlated with gas surface density and metallicity. 
The correlation is consistent with the mean extinction depending on dust density 
where the dust to gas mass ratio scales with the metallicity. This trend is expected if 
H II regions tend to be located near the mid-plane of a gas disk and emerge from their 
parent molecular clouds soon after birth. In environments with gas densities below a 
few hundred M0/pc^, star formation rates estimated from integrated line fluxes and 
mean extinctions are likely to be fairly accurate. 



Subject headings: 



1. Introduction 

It has long been suspected that H II regions tend to exhibit higher reddenings in higher 
metallicity environments. For example, in MlOl and M51 ( ^cowen, Dufour, fc Hester 1992 ; 



Kennicutt & Garnett 19961; Scowen 1992| ) a trend between H II region extinction and galact 



ic 



radius was observed and interpreted to be caused by the radial metallicity gradient in these 
galaxies. However, most studies have found a large scatter between the extinction and metallicity 



^Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science 
Institute which is operated by the Association of University for Research in Astronomy, Inc. (AURA), under NASA 
contract NAS5-26555. 



- 2 - 



in individual H II regions ( McCall, Rybski fc Shields 1985| ; Zaritsky, Kennicutt & Huchra 1994; 
Scowen et al. 1992; Kennicutt & Garnett 1996; Bresolin et al. 1999D . Even though the metallicity 
is fairly uniform across the Large Magellanic Cloud (LMC), high values of the Balmer decrement 
are measured near 30 Doradus ( Heydari-Mayaleri fc Testor 1985| ; Caplan & Deharveng 1985) 
which is located in a region of high gas density (Cohen et al. 1988; Luks fc Rohlfs 1992| ). ^cowen 



et al. 1992| discovered a few H II regions in MlOl with anomalously high extinctions at large 
galactic radii, which they attributed to higher local gas densities in nearby spiral arms. Despite 
the large scatter intrinsic to ensembles of H II regions, previous studies have found that the 



mean extinction of ensembles of H II regions varied between different galaxy regions. Martin fc 



Friedli 1999 attributed the higher mean extinctions in H II regions in barred galaxies compared to 
non-barred galaxies to be caused by the higher gas densities associated with bar induced shocks. 



Scowen et al. 1992 and Scowen 1992 found that mean H II region redennings were higher at large 



galactic radii than at smaller radii in both MlOl and M51. 

If there is a general trend for an ensemble of H II regions to be more highly reddened in 
higher gas density or metal rich environments, then the star formation rate estimated from a sum 
of recombination line (such as Ha; e.g., Kennicutt 199^ ) fluxes would be underestimated in these 
environments. One way to quantify this effect is to measure extinctions towards ensembles of H II 
regions in a variety of different galactic environments. If mean H II region extinctions depend upon 
gas density and metallicity then it might be possible to formulate an empirical correction for the 
flux. This would allow star formation rates to be measured from one recombination line in more 
heavily reddened regions. If such a trend exists, star formation rates could be more accurately 
estimated for galaxies at high redshift and observed line ratios would also yield constraints on the 
galaxy metallicity and gas surface density. 

By observing in the the near-infrared it is possible to probe into more heavily obscured regions 
because the near-infrared is less affected by dust extinction than optical wavelengths (for examples 
of previous work on infrared recombination lines see Calzetti et al. 1996 and Beck, BeckwithJ 
fc Gatley 198^ ). The hydrogen Paa recombination line at 1.875/im is bright in the vicinity of 
hot newly formed stars. It is difficult to observe in nearby galaxies from the ground because of 
strong atmospheric absorption from water molecules. It is expected to be about one eighth as 
bright as Ha at temperatures and densities typical of H II regions, but since it is three times 
the wavelength of Ha it should be less affected by extinction. For an H II region seen through a 
moderate foreground dust screen with Ay ~ 3, Paa should be as bright as Ha. Previous studies 
of Paa imaging in galaxies include Boker et al. 199!^ and Marconi et al. 2000| . 



In this paper we measure Ha and Paa fluxes for H II regions in the central regions of five 
galaxies. We find that the scatter in Ay in each galaxy region is less than the differences in the 
mean values. Observations of distant galaxies are restricted towards observing line emission over a 
large region in the galaxy. Because of this we estimate a correction for the Ha flux which is based 
upon summing over the extinction corrected fluxes of individual H II regions. This procedure is 
one way of averaging over the scatter in extinctions which has been observed in previous studies 
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of H II regions. We then investigate the dependence of this mean extinction on local gas density 
and metallicity. 



2. Observations 

We searched the HST archive for galaxies which have high quality narrow band images in Paa 
and reveal a number of H II regions. The archived images are listed with their filters in Table 1. 
NGC 2903 and NGC 1512 have narrow band images taken by HST in both Paa and Ha, revealing 
a large number of H II regions. M51 also has both Paa and Ha images observed by HST showing 
numerous H II regions; however, these images overlap only in a small region about 20" to the 
north and 20" to the east of the nucleus. NGC 4449 and NGC 6946 have high quality Paa images 
containing numerous H II regions but lack HST narrow band imaging in Ha. For comparison we 
used Ha fluxes for NGC 4449 measured by ^cowen 1992 and a calibrated ground-based Ha image 



for NGC 6946 provided by Rob Kennicutt (private communication). NGC 2903, NGC 4449 and 
NGC 6946 were observed as part of the nearby galaxies snap shot survey ( iBoker et al. 19991) . 



The NICMOS images were reduced with the nicred data reduction software ( [VIcLeod 1997] ) 
using on orbit darks and flats. The narrow band images required a pedestal correction which we 
assumed was an additional constant in each quadrant of the chip (for more discussion on this see 
Boker et al.T999| ). WFPC2 images were combined using the IRAF/STSDAS routine combine 



which removed most of the cosmic ray hits. We derived Ha and Paa line emission images by 
subtracting scaled continuum images from the narrow band images containing the lines. 

WFPC2 images were calibrated using zeropoints listed in the HST Data Handbook. 
Specifically we used 3.062 x 10~^^ and 2.941 x 10~^^ erg cm~^ s~^/(DN s~^) for the continuum 
subtracted F656N and F658N narrow band images, respectively in the PC camera of WFPC2. 
The NICMOS continuum subtracted narrow band images in the F187N filter were calibrated using 
6.69 X 10^16 and 8.27 x 10"^^ erg cm ^ s ""^/(DN s ^) for Camera 2 and Camera 3 respectively 
for the Paa line. The NICMOS broad-band images in the F160W filter were calibrated using 
2.190 X 10"^ and 2.776 x 10~^ Jy/(DN s"^) for Camera 2 and Camera 3 respectively. The 
NICMOS flux calibration is 10.7% lower than that used by [Boker et al. 1999 and is based on 



measurements of the standard stars P330-E and P172-D during the Servicing Mission Observatory 
Verification program and subsequent observations (M. Rieke 1999, private communication). The 
filter widths listed in the NICMOS and WFPC2 instrument handbooks were used to derive the 
calibration corrections. The location of the lines was near the region of maximum transmission in 
the filters in all cases so no correction for filter transmission is required. 

To measure the flux of emission regions we used aperture photometry. We identify individual 
regions from peaks in the Paa line emission images. Aperture sizes were chosen to match the 
resolution of the Ha images or that of existing Ha photometry (for NGC 4449 to match that 
of ^cowen 199^ ). For NGC 2903, NGC 1512, and M 51, because we had HST imaging in both 
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lines we were able to use very small apertures. In all cases identical aperture positions and sizes 
were used on both line emission images. Fluxes and identified regions are listed in Tables 2-6. 
Positions are given with respect to the nucleus as defined by the centroid at the brightest spot 
in the NICMOS/F160W images. The positions of nuclei are listed in Table 2-6. Astrometry was 
done with the IRAF/STSDAS routine xy2rd from the F160W images and is that determined by 
the HST pointing. 

We use the Ha/Paa line ratios to estimate the extinction to each region. We assume intrinsic 
flux ratios of Fho/Fhis = 2.86, Fna/Fpaa = 8.46 for Case B recombination at a temperature of 
lO^K and a density 100 cm^'^ (line ratios are given in psterbrock 1989 ) and a Galactic extinction 
law (pVIathis 19901) . For NGC 2903 and M51 the [NII]6584l line was not contained in the F656N 
filter, however it is contained in the F658N image of NGC 1512. We correct for the presence of 
this line in NGC 1512 and NGC 6946 assuming a Ha/(Ha + [Nil]) hue ratio of 0.75 (following 
the procedure of Kennicutt 1983). We corrected for Galactic redenning using values from ^chlegel 
et al. 199^ . 



Fluxes and extinctions derived from Paa/Ha flux ratios are listed in Tables 2-6. The 
variations in the intrinsic line ratios caused by temperature differences are unlikely to significantly 
affect the measured extinctions; (lowering the electron temperature by 5000K results in a reduction 
of of 15% in the Ha/Paa line ratio). The errors in the measured line ratios are dominated by 
errors in the sky value caused by color variations which affect our continuum subtraction or the 
uneven sky in the Paa images caused by the pedestal affect. The photometric uncertainties are 
roughly estimated in the Table captions. There is an additional error in the absolute calibration 
of ~ ±20% which corresponds to an error of ~ ±0.2 in the derived Ay- Relative photometric 
error between individual H II regions in a galaxy are ~ ±0.1 mag in Ay for the faintest regions. 
The scatter, ~ ±0.5, in Ay in each galaxy is therefore likely to be intrinsic and not caused by 
photometric uncertainty. 

When possible, we have compared our measured extinction values with those reported in 



previous studies. Two of the H II regions in NGC 6946 were also studied by Byman et al. 200C 
with radio continuum and Ha observations. For the 3rd and 4th H II region listed in Table 5 we 
find Ay ~ 3.2, while Hyman et al. 200C report Ay ~ 2.5 and 3.4, in rough agreement with our 
estimates. 



Although Scowen 1992 measured a high extinction. Ay ~ 2, for many regions in NGC 4449 
from his Ha and H/? narrow band imaging, the Paa/Ha line ratios that we compute from our 
Paa photometry and his Ha fluxes suggest that there is little internal extinction in this galaxy. 



However, our line ratios are consistent with the spectroscopic measurements of Kobulnicky, 



Kennicutt, & Pizagno 1999 who find Ha/H/3 ~ 2.86 which means there is little extinction. 
We suspect that there may have been an error in the calibration of the H/5 image by pcowen 



1992. Ha fluxes listed by Fuentes-Masip, Castaneda, & Munoz-Tunon 2000 are consistent with 



those measured by Scowen 1992 so the Ha fluxes measured by Scowen 1992 are sufficiently 
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well-calibrated. It is a concern that extinctions measured from the Balmer decrement might be 
biased by H/3 absorption present in the stellar component (e.g., Dufour et al. 1980 ). However, 
after correcting for H/3 absorption from the stellar component, pcowen 1992| found no strong 
changes in his correlation plots. 



2.1. Comparing Ha and Paa Morphology 

In Figures 1-3 we show the HST images for NGC 1512, NGC 2903, M 51 in both Paa and 
Ha lines. We found little correspondence between star clusters observed in the broad band images 
and the line emission maps. This makes it difficult to constrain the ages of the cluster based on 
Paa equivalent widths. A number of star clusters were observed which were outside the region 
of line emission emission and so could correspond to regions where radiation pressure, ionization 
fronts, winds and supernovae have cleared the clusters of gas. In NGC 1512 and NGC 4449 the 
luminosities in Ha of the faintest regions are ~ 10^^ ergs/s and so similar in luminosity to the 
Orion nebula ( Kennicutt 1984] ) which contains one O star. Since H II regions cannot be much 



fainter than this, we suspect that there are no deeply embedded H II regions in these two galaxies. 
The faintest regions detected in NGC 2903 and M51 are an order of magnitude brighter and so 
could contain a few dozen O stars. 

In NGC 1512, NGC 2903 and M51, there are few qualitative differences between the respective 
line emission images. In other words, almost all H II regions above the detection limit in Paa were 
seen in both images. One region stands out as an exception to the above statement. This region 
in NGC 2903 is 0'.'71 east and 4'.'17 north of the nucleus (see Figure 1) and is significantly brighter 
in Paa than in Ha. From the line ratio we estimate Ay ~ 4, however we did not see a peak in the 
Ha line image; instead the line emission is diffuse and may not be emitted from the same spatial 
region. Previous studies have measured large extinctions in the nucleus of this galaxy (Ay ~ 15; 
Lebofsky fc Rieke 1979 ) as well as high gas densities ( packson et al. 1991 ) so this region may be 



an example of an obscured region. 

Total Paa and Ha fluxes across our NGC 2903 images are 4.5 x 10-^3 and 5.3 x lO-^^ 
ergs cm~^s~^ respectively. We can compare these values to the Bra and Br7 fluxes measured 



by Beck, Beckwith, & Gatley 1984 in an 8" diameter aperture. We would predict Bra and 
Br7 fluxes of 1.53 x 10^^^, and 5.30 x 10^^^ ergs cm~^s~^ respectively, using an Ay of 3.06 
consistent with the Paa and Ha line ratio across this aperture. Our estimated Bra and Br7 fluxes 



are consistent with those measured by [Beck, Beckwith, fc Gatley 1984 , 1.4 it 0.18 x 10 and 
3.5 it 1.2 X 10~^^ respectively. The extinction we estimate Ay ~ 3 — 4 is therefore consistent with 



line fluxes and ratios measured by Beck, Beckwith, &: Gatley 1984, but substantially lower than 



that. Ay ~ 15, measured from broad and narrow band imaging near the 10/im silicate absorption 
feature ( [Lebofsky fc Rieke 1979|) . The region identified above is not bright enough (it is only 



l/20th the total) in Paa to dominate the total Paa emission in the 8" nuclear region. Since the 
line ratios agree with Ay ~ 3 — 4, it is also unlikely that this region dominates the Bra emission. 
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The total 10/im flux roughly corresponds to that expected if a moderate fraction (e.g. 1/2) of the 
total bolometric luminosity in the UV (needed to account for the line emission) is re-radiated by 
dust ( Peck, Beckwith, fc Gatley 1984|) . The 10/xm emission is most likely to be emitted from the 
extended region we see in the line emission maps and not from a single deeply embedded region. 
The most likely explanation for the discrepancy between the extinctions we measure and that 
estimated at 10/xm is that the 10/im spectrum contains structure from by PAH emission which 
caused Lebofsky &: Rieke 1979 to overestimate the depth of the silicate absorption feature (see 
for example ^iebenmorgen, Kriigel, fc Zota 1999]) . The extinction we estimate near the nucleus of 



NGC 6946 {Ay ~ 4.6) is also significantly lower than that estimated by Lebofsky &: Rieke 1979 
{Ay ~ 50). 



3. A comparison between gas surface density, metallicity and extinction 

For a distant galaxy we would like to to derive star formation rates based on the emission line 
flux. We therefore desire a way to estimate the total emission in this line that would be emitted in 
the absence of extinction. In the comparison between Ha and Paa images we found no evidence 
for a population bright HII regions that are extremely red compared to the average. The scatter 
in Ay in each galaxy region (approximately ±0.5) is less than the differences in the mean values 
which range from a mean Ay ~ 3.3 to 0. Because the mean values span a large range, we are 
prompted to look for correlations between this value and environmental properties. 

To compare extinctions in different galactic environments we must first compile information 
from the literature on the metallicities and gas densities in the various regions of these galaxies. 
These are listed in Table 7 with references to the literature from which they were taken. The 
gas densities are primarily taken from single dish CO measurements with similar aperture sizes. 
When necessary we have corrected published values of CO so that we used the same conversion 
factor of N(H2) = 2.8 x lO'^^cm-'^Ico (K km s'^) (Scoville et al. 1987). The metallicities were 
primarily taken from Zaritsky et al. (1994). We estimate that the uncertainties in the gas density 
(approximately a factor of 2) and metallicity measurements (approximately 0.3 dex) are smaller 
than the variations between the measurements in different galaxies and galaxy regions. 

We did not find a measured gas mass in NGC 1512 so we estimate the gas surface density 
based on the color map made from the broad band V (0.5;um ) and H band (1.6//m) images. We 
observed a change in color of A(V-H) 0.4 between the circumnuclear ring and the region interior 
and exterior to it. We assume that the dust is a foreground screen for half of the light (for the 
bulge stars that are located behind the gas). This is equivalent to assuming that the scale length 
of the gas is smaller than that of the stars. This change in color then corresponds to Ay ~ 1.3 and 
a gas surface density in hydrogen of T^h ~ 19M0 pc~^ (based on the Ay to N(H) conversion and 



extinction law of Mathis 1990 ). 



We did not find any optical spectroscopy for NGC 1512 in the literature. However, we can use 
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the correlation between absolute magnitude and central metallicity found by parnett et al. 1997 
to estimate the metallicity in its circumnuclear ring. This galaxy has an absolute magnitude of 
Mb ~ —18.9 (assuming a distance of 9.5 Mpc from [Tully 1988 with a Hubble constant of 75 
km Mpc~^). The galaxy has an inclination corrected circular velocity of Vc ~ 173 km/s. From 
Figure 4 of Garnett et al. 1997] we estimate that NGC 1512 should have a central metallicity that 
is approximately solar. 

To increase the size of our comparison sample we also compile the properties of H II regions 
as a function of radius in MlOl and M51. H II region extinctions and luminosities in MlOl were 



taken from pcowen et al. 1992| and are restricted to the 200 brightest H II regions. To describe 
the azimuthal gas density in MlOl we summed the azimuthal molecular and atomic hydrogen gas 



mass densities reported by Kenney, Scoville & Wilson 1991. The metallicity gradient was taken 



from Kennicutt & Garnett 1996. H II region extinctions and luminosities in M51 were taken from 



Scowen 1992 and are restricted to be from the brightest H II regions defined as having an observed 
flux in Hq greater than 10"^"^ erg cm~^ s~^. Molecular and atomic hydrogen gas densities in M51 
were taken from pydbeck, Hjalmarson &: Rydbeck 1985| and Scoville & Young 1983. The oxygen 



abundance gradient in M51 were taken from Zaritsky, Kennicutt & Huchra 1994 



3.1. Weighting 

We approach the problem of correcting for extinction by considering a sum over individual H II 
regions and by deriving a correction that can be multiplied by the total observed flux in the line 
to obtain the true one. This procedure is one way of averaging over a large scatter in extinctions 
and allows us to look for trends between the mean extinction and environmental properties. The 
correction we can describe as a luminosity weighted mean extinction. To derive this correction 
we sum over the corrected fluxes of each region, so it involves a luminosity weighted sum. This 
makes the estimate more robust than using the mean extinction computed from the distribution 
of extinctions measured towards each region. The mean computed from the distribution would 
be more dependent on the faintest HII regions and so noisier and less robust than the luminosity 
weighted mean extinctions computed as in equation (3) below. 

We denote the flux in Ha from a region i by FHa,i- We denote the sum of observed fluxes 
over a larger area as Fna = Si FHa,i- We denote the true or corrected emission as 

a,i,corr- (1) 

i 

We aim to estimate a factor ana such that we can multiply the observed total flux in a large 
region to obtain the true or corrected emission. 

FHa,corr = aHaFua- (2) 

We correct the fluxes from each region using the extinctions derived from the line ratios and then 
sum the regions to estimate FHa,corr- At the same time we sum the uncorrected fluxes to get Fna- 
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From the ratio of Fjja and FHa,corr we estimate the correction factor ana that we need to obtain 
the true or corrected emission from the observed total. This is converted to a mean extinction at 
the wavelength of Ha by 

A{lQ=-2.5logaHa. (3) 

A similar mean extinction can also be calculated for Paa. This can be compared to the extinction 
estimated from the ratio of the total Paa and Hq fluxes which we refer to as Atotai{Ha). 

We calculate the luminosity weighted mean extinction (defined above) by summing over 
the Paa and Ha fluxes listed in Tables 2-6 for the nuclear regions of NGC 1512, NGC 2903, 
NGC 4449, M51 and NGC 6946. These mean extinctions are listed in Table 7. We also perform 
these sums using Ha and H/3 fluxes as a function of radius in MlOl and M51 based on the Balmer 
decrement. The radial width over which we summed the properties of H II regions was 1 kpc for 
M51 and 1.5 kpc for MlOl. For comparison to an additional low metallicity system we also derive 
a similar estimate for the LMC based on the Balmer decrement. Ha and H/3 fluxes for LMC H II 



regions were taken from Caplan fc Deharveng 1985 and we correct for Galactic extinction using a 



value of Ab = 0.324 (|Schlegel et al. 199g| ). 



In Table 7 we can compare A{Ha) with Atotai{Ha). A{Ha) allows us to correct the total 
observed Ha flux for redenning based on the reddenings to each H II region. We see a tendency 
for A(Ha) to be somewhat higher than that estimated from the line ratio of the integrated 
fluxes {Atotai{Ha)). The differences are small because the dispersion in the extinctions between 
individual H II regions in each galaxy is low. A star formation rate estimated from the total Ha 
line flux and corrected with an extinction based on the Ha/Paa line ratio is likely to be quite 
accurate. 

We plot our mean extinctions (defined in Eqn 3) versus gas density in Figure 4a. As expected 
there is a correlation between mean extinction and local gas density. However some of the 
correlation may be due to the dependence of the CO to H2 conversion factor on metallicity. To 
take this dependence into account we correct the H2 by a factor that depends on the metallicity 
and display the corrected gas density in Figure 4b. Our correction assumes an increase in the 
conversion factor by 4.6 when the oxygen abundance drops by a factor of 10 ( [Wilson 1995 ). 



Extending this relation towards metallicities larger than solar is reasonable since galaxies with 



metallicities greater than solar are expected to have lower conversion factors (e.g., Nakai & Kuno 



1995). The correlation between gas density and mean extinction is not significantly changed. 



It is difficult to separate between the role of gas density and that of metallicity since the 
densest regions also correspond to the most metal rich. However as we see in Figure 4c, metallicity 
alone is not sufficient to account for the large variation in our mean extinctions. 

The dust density is the quantity most likely to affect the reddening so we desire a comparison 
between our extinction values and the dust surface density. We can estimate the dust surface 
density from the gas density using a metallicity dependent gas to dust mass ratio. The gas to 
dust ratio may correlate better with the carbon abundance than with the oxygen abundance 
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(Bohlin, Savage, & Drake 1978; Martin, Maurice, fc Lequeux 19891 ) . To convert between 



oxygen abundance and carbon abundance we use the relation measured by Garnett et al. 1995 
[log(C/0) « 0.54 log(0/H) +1.6]. We assume that at solar metallicity the gas to dust mass 
ratio is 120. To estimate the dust surface density we multiply the solar dust to gas ratio by the 
gas density and the ratio between the carbon abundance predicted from the oxygen abundance 
and that at solar. The relation between mean extinction and estimated dust density is shown 
in Figure 4d. The dashed line plotted in Figure 4d corresponds to that predicted for a uniform 
medium with the H II regions located in the middle of the dust distribution. We find that the 
correlations we have seen between mean extinction, gas density and metallicity are consistent with 
an approximately linear relation between the log of the dust density and mean extinction. 

For the high dust densities, the mean extinctions expected are somewhat higher than those 
observed. In Figure 4d we plotted the extinction predicted if the H II regions were in the mid plane 
of a uniform disk. These deviations would be predicted from a model which takes into account 
the patchiness of the ISM (e.g., Calzetti et al. 1996| ). To illustrate this we have also plotted in 
Figure 4d a patchy dust model for a two component dust model with an extinction ratio of four 
and an area filling factor of 0.5. However such a model would result in a larger scatter (greater 
than ±0.75) in the extinctions towards individual regions than we measured. A more realistic 
patchy dust model may match the extinction and mean gas density. There could be nearby dense 
molecular clouds lacking H II regions which would have contributed to the mean gas density in the 
region but not to the H II region extinctions. 



4. Summary and Discussion 

In this paper we have presented a comparative study between Paa and Ha in the central 
regions of five galaxies. For three galaxies we have HST images in both lines and the morphologies 
of the respective line emission maps are remarkably similar. The dispersion in extinctions is about 
a magnitude for the H II regions in each galaxy and is smaller than the differences between the 
mean extinctions. We find no evidence for a population of bright heavily embedded H II regions. 

Only one bright region, in NGC 2903, was detected clearly in the Paa image and was faint or 
not seen in the Ha image. However this region does not dominate the Paa flux of the nuclear star 
forming region and is also unlikely to dominate the Bra emission. This suggests that an estimate 
of the star formation rate based on a mean extinction and either the Ha or Paa flux is likely to be 
fairly accurate. We suspect that some previous studies of the silicate absorption feature at 10/um 
may have overestimated the extinction towards the nuclei of NGC 2903 and NGC 6946 because of 
contamination by PAH features. 

We have compared mean extinctions as a function of gas density and metallicity and find 
correlations between these quantities. The correlations between mean extinction, gas density 
and metallicity are consistent with an approximately linear relation between the log of the dust 
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density and mean extinction. This trend is expected if the H II regions tend to be located near the 
mid-plane of a gas disk and emerge from their parent molecular clouds soon after birth. Radiation 
pressure and ionization fronts are expected to clear a young star forming region of dense gas on a 
timescale of 10^ years ( [Ferland 1998| ) whereas H II regions should exist for a few times 10^ years. 
It would be unlikely to detect an H II region in a region with extinction significantly higher than 
the mean predicted from the local gas density. This may provide a partial explanation for the low 
dispersion in extinctions that we measured in each galaxy. 

We expect that forthcoming high resolution CO and HI surveys will allow a more careful 
comparison between the statistical properties of H II regions and the gas distribution. It may be 
possible to study correlations between the patchy distribution of the ISM, the tendency for star 
forming regions to evacuate themselves of dense gas, and large scale gaseous structures cause by 
bars and spiral arms. The relationships presented here may allow more accurate estimates of 
star formation rates in galaxy regions which contain gas densities up to a few hundred Mq pc~^. 
Above this level it is likely that a significant fraction of the line emission is absorbed even in the 
near-infrared and mid or far infrared measurements would be required (e.g., Kennicutt 19981 ). 
Further work on the trends discussed here may make it possible to more accurately measure star 
formation rates and constrain gas densities in distant galaxies. 
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Fig. 1. — HST images of the circumnuclear ring of NGC 2903. The Ha image was constructed 
WFPC2 images. The Paa was observed with NICMOS Camera 3 and constructed from F187N and 
F160W images. The V-H color map was constructed from the F555W and F160W images. Colors 
from V-H = 2.0 (white) to 3.5 (black) are shown. Axes are given in arcsec from the nucleus. There 
is a bright source at 0'.'71 east, and 4'.'17 north from the nucleus evident in the Pact image but not 
the Ha image. 
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Fig. 2. — HST images of the circumnuclear ring of NGC 1512. The Ha image was constructed 
WFPC2 images. The Paa was observed with NICMOS Camera 2 and constructed from F187N and 
F160W images. The V-H color map was constructed from the F547M and F160W images. Colors 
from V-H = 2.5 (white) to 3.3 (black) are shown. Axes are given in arcsec from the nucleus. We 
found no evidence for bright embedded H II regions. 
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Fig. 3. — HST images of a region about 20" to the north and 20" to the east of the nucleus in 
M51. The Ha image was constructed from WFPC2 images. The Paa was observed with NICMOS 

Camera 3 and constructed from F187N and F190N images. The V-H color map was constructed 
from the F547M and F160W images. Colors from V-H = 2.5 (white) to 3.3 (black) are shown. 
Axes are given in arcsec from the brightest region listed in Table 6. See Table 6 for the coordinates 
of this region. 
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Fig. 4. — a) The relation between the luminosity weighted extinction and the azimuthal hydrogen 
gas density. NGC 1512, NGC 2903, NGC 4449, NGC6946, the LMC, and the north east position 
of M51 are plotted by name. The MlOl and M51 azimuthal averages are plotted as solid triangles 
and open squares, respectively, where each point represents an average at a different radius. The 
dashed line shows the linear relation between Ay and half of the gas column depth (assuming the 
relation between hydrogen column depth and extinction from Mathis 1990). Wc use half of the 
gas column, assuming that the H II regions are embedded in the gas disk. The gas densities have 
been estimated with a CO conversion factor of N(H2) = 2.8 x 10^*^cm""^/co (K km s~^) (Scoville 
et al. 1987). b) Similar to a) however we have corrected the hydrogen gas densities by using a 
CO conversion factor that depends on metallicity (see text for details), c) Dependence of the 
luminosity weighted extinction with metallicity only. We see this correlation in part because the 
regions with higher gas density also correspond to regions of higher metallicity. The dotted line 
shows the extinction predicted from a gas surface density of IOMqpc"^ and a dust to gas fraction 
that depends on the carbon abundance. The metallicity alone is not sufficient to account for the 
variations in extinctions, d) Dependence of the luminosity weighted extinction (defined in Eqn 3) 
with estimated dust surface density. The dust surface density was estimated from the corrected gas 
densities (shown in b) by scaling with the carbon abundance (see text for details). The dashed line 
shows the linear relation between Ay and half of the dust column depth. The dot dashed line is a 
2 component patchy dust model with an area filling factor of 50% and an opacity ratio of 4. This 
model illustrates that a patchy medium would have a lower mean extinction than that predicted 
directly from the dust density. 




10 -10 10 -10 



Fig. 1. — HST images of the circumnuclear ring of NGC 2903. The Ha image was constructed 
WFPC2 images. The Paa was observed with NICMOS Camera 3 and constructed from F187N and 
F160W images. The V-H color map was constructed from the F555W and F160W images. Colors 
from V-H = 2.0 (white) to 3.5 (black) are shown. Axes are given in arcsec from the nucleus. There 
is a bright source at 0."71 east, and 4." 17 north from the nucleus evident in the Paa image but not 
the Ho; image. 




Fig. 2. — HST images of the circumnuclear ring of NGC 1512. The Ha image was constructed 
WFPC2 images. The Paa was observed with NICMOS Camera 2 and constructed from F187N and 
F160W images. The V-H color map was constructed from the F547M and F160W images. Colors 
from V-H = 2.5 (white) to 3.3 (black) are shown. Axes are given in arcsec from the nucleus. We 
found no evidence for bright embedded Hll regions. 




Fig. 3. — HST images of a region about 20" to the north and 20" to the east of the nucleus in 
M51. The Ha image was constructed from WFPC2 images. The Paa was observed with NICMOS 
Camera 3 and constructed from F187N and F190N images. The V-H color map was constructed 
from the F547M and F160W images. Colors from V-H = 2.5 (white) to 3.3 (black) are shown. 
Axes are given in arcsec from the brightest region listed in Table 6. See Table 6 for the coordinates 
of this region. 



a 



1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 ,1 1 1 1 1 1 1 




/ N6946 




q □ 








/ 

, N2903 




/ M51NE 






- 




N1512 ° 


' D 


□ ' 




□ / 








A/ 




^ □ 








LMC 








N4449-''V ^ 

" , 1 , , , , 1 , , , , 1 , 


, , , 1 , , , , 1 , , " 


0.5 1 1.5 


2 2.5 



log(S^) Msol pc 



X 2 



1-771- 

f 

ts(fe946 



N2963 
/M51NE 



N1512 /□ 



LMC 



0.5 1 1.5 2 2.5 

log(Epj) Corrected Msol pc"'^ 



N6946 

□ □ 



N2903 
M51NE 



-tMCT 

N444^ A * 



_L 



8.5 9 9.5 

12 + log(0/H) 



X 2 



M2903 

M51NE_ 



-1.5 -1 -0.5 
log(^duJ Msol pc-^ 



Fig. 4. — a) The relation between the himinosity weighted extinction and the azimuthal hydrogen 
gas density. NGC 1512, NGC 2903, NGC 4449, NGC6946, the LMC, and the north east position 
of M51 are plotted by name. The MlOl and M51 azimuthal averages are plotted as solid triangles 
and open squares, respectively, where each point represents an average at a different radius. The 
dashed line shows the linear relation between Ay and half of the gas column depth (assuming the 
relation between hydrogen column depth and extinction from Mathis 1990). We use half of the 
gas column, assuming that the HII regions are embedded in the gas disk. The gas densities have 
been estimated with a CO conversion factor of N(H2) = 2.8 x 10^''cm~^/c'o (K km s~^) (Scoville 
et al. 1987). b) Similar to a) however we have corrected the hydrogen gas densities by using a 
CO conversion factor that depends on metallicity (see text for details), c) Dependence of the 
luminosity weighted extinction with metallicity only. We see this correlation in part because the 
regions with higher gas density also correspond to regions of higher metallicity. The dotted line 
shows the extinction predicted from a gas surface density of IOMqpc^^ and a dust to gas fraction 
that depends on the carbon abundance. The metallicity alone is not sufficient to account for the 
variations in extinctions, d) Dependence of the luminosity weighted extinction (defined in Eqn 3) 
with estimated dust surface density. The dust surface density was estimated from the corrected gas 
densities (shown in b) by scaling with the carbon abundance (see text for details). The dashed line 
shows the linear relation between Ay and half of the dust column depth. The dot dashed line is a 
2 component patchy dust model with an area filling factor of 50% and an opacity ratio of 4. This 
model illustrates that a patchy medium would have a lower mean extinction than that predicted 
directly from the dust density. 



Table 1. HST archival images 



galaxy Paa continuum Ha continuum 



NGC 1512 F187N/NIC2 F160W/NIC2 F658N F547M 

NGC 2903 F187N/NIC3 F160W/NIC3 F656N F555W 

M51 F187N/NIC3 F190N/NIC3 F656N F547M 

NGC 4449 F187N/NIC3 F160W/NIC3 

NGC 6946 F187N/NIC3 F190N/NIC3 



Note. — Line and continuum filters on HST used to make the 
Paa and Ha line emission maps shown in Figures 1-3 and used to 
measure H II region intensities and line ratios listed in Tables 2-6. 
The near-infrared images were observed with NICMOS and the filter 
names are followed by the NICMOS camera used for the observation. 
The pixel size for the NICMOS camera 2 (NIC2) is ~ 0'.'076, and for 
the NICMOS camera 3 (NIC3) is ~ 0'.'204. 



Table 2. H II regions in NGC 1512 
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Note. — We did aperture photometry in 0^'75 radius apertures. Offsets 

are given with respect to the position of the nucleus; RA(2000) = 4:03:54.131, 
DEC(2000)= -43:20:56.54, as measured from the HST pointing. Columns: (1) 
Offset in RA from the nucleus. A positive offset is east of the nucleus; (2) Offsets 
in DEC from the nucleus. A positive offset is north of the nucleus; (3) Measured 
flux in Ha. Uncertainty in the sky value causes an error of ~ ±5% in the faintest 
regions; (4) Ratio of flux in Paa to Ha. Uncertainty in the sky value causes 
an error of '~ ±13% in the measurement of the Paa flux in the faintest regions; 
(5) Ay computed from the ratio. The distribution of extinctions is centered at 
Ay = 1.8 with a standard deviation of 0.7; (6) Luminosity in Ha which has been 
corrected for extinction using Ay given in column (5). We assume a distance of 
9.5 Mpc (Tully 1988; Hq = 75 km Mpc^^). 



Table 3. H II regions in NGC 2903 
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Note. — Wc did aperture photometry in 0'.'469 radius apertures. Offsets 
are given with respect to the position of the nucleus; RA(2000) = 9:32:09.98, 
DEC(2000)= -^21:30:0.8, as measured from the HST pointing. Columns: (1) 
Offset in RA from the nucleus. A positive offset is east of the nucleus; (2) Offset 
in DEC from the nucleus. A positive offset is north of the nucleus; (3) Measured 
flux in Ha. Uncertainty in the sky value causes an error of ~ ±2% in the faintest 
regions; (4) Ratio of flux in Paa to Ha. Uncertainty in the sky value causes an 
error of ~ ±4% in the measurement of the Paa flux in the faintest regions; (5) Ay 
computed from the ratio. The distribution of extinctions is centered at Ay = 3.3 
with a standard deviation of 0.7; (6) Luminosity in Ha which has been corrected 
for extinction using Ay given in column (5). We assume a distance of 6.3 Mpc 
(TuUy 1988; Hq = 75 km s'^ Mpc'^). 



Table 4. H II regions in NGC 4449 
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Note. — We did aperture photometry in 4" radius apertures to match the apertures 
used by Scowen (1992) to measure Ha fluxes. Columns: (1) RA of emission peak in RA as 
we measure in the Paa image; (2) DEC of emission peak as measured in the Paa image; 
(3) Measured flux in Paa. Uncertainty in the sky value causes an error of ~ ±11% in the 
measurement of the Paa flux in the faintest regions; (4) Ratio computed by using our Paa 
flux and the Ha fluxe measured by Scowen (1992) in the same region; (5) Ay computed 
from the ratio; (6) Luminosity in Ha which has been corrected for extinction using given 
Ay given in column (5). We assume a distance of 3.05 Mpc (Stewart et al. 2000). 



Table 5. H II regions in NGC 6946 
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Note. — We did aperture photometry in 2'.'7 radius apertures. Offset are given 
with respect to the position of the nucleus; RA(2000) = 20:34:52.33, DEC(2000)= 
-1-60:09:14.21, as measured from the HST pointing. Columns: (1) Offset in RA 
from the nucleus. A positive offset is east of the nucleus; (2) Offset in DEC for 
the nucleus. A positive offset is north of the nucleus; (3) Measured flux in Paa. 
Uncertainty in the sky value causes an error of ~ ±11% in the measurement of 
the Paa flux in the faintest region; (4) Ratio of flux in Paa to Ha. Uncertainty in 
the sky value causes an error of ~ ±15% in the measurement of the Ha flux in the 
faintest region; (5) Ay computed from the ratio. The distribution of extinctions 
is centered at Ay = 2.8 with a standard deviation of 1.2; (6) Luminosity in Ha 
which has been corrected for extinction using Ay given in column (5). We assume 
a distance of 5.5 Mpc (Tully 1988; Hq = 75 km s"^ Mpc"^). 



Table 6. H II regions in M51 
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Note. — We did aperture photometry in 0'.'41 radius apertures. Offsets are 
given with respect to the first listed position; RA(2000) = 13:29:54.00, DEC(2000) 
= 47:12:11.86. Columns : (1) Offset in RA from the first listed position. A 
positive offset is east of the first position; (2) Offset in DEC from the first position. 
A positive offset is north of the first detected position nucleus; (3) Measured flux 
in Ha. Uncertainty in the sky value causes an error of ~ ±4% in the measurement 
of the Ha flux in the faintest region; (4) Ratio of flux in Paa to Ha. Uncertainty 
in the sky value causes an error of ~ ±8% in the measurement of the Ha flux in 
the faintest region; (5) Ay computed from the ratio; (6) Luminosity in Ha which 
has been corrected for extinction using Ay given in column (5). The assumed 
distance is 9.6 Mpc (Sandage & Tammann 1975). 



Table 7. Integrated properties of H II regions, local gas density and metalicity 
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Note. — Values for the mean extinctions A{Ha) and A{Paa) (defined in Equation 
3) were derived from the Paa fluxes and Paa/Ha flux ratios listed in Tables 2-6 for all 
galaxies except the LMC. This can be compared to the extinction estimated from the ratio 
of the total Paa and Ha fluxes which we refer to as Atotal{Ha). The value for A{Ha) 
for the LMC was estimated from the fluxes and Balmer decrements listed by Caplan & 
Deharveng (1985). Metalicities are estimated for the region covered by our images. For 
NGC 1512 the central metalicity was estimated based on its luminosity and a gas density 
corresponding to the extinction we observe in our color map (see text). Central metalicities 
for NGC 2903, M51, and NGC 6494 were estimated from gradients and values listed by 
Zaritsky et al. (1994). The metalicity of NGC 4449 is that measured by Heckman et al. (1988). 
The gas density for the nucleus of NGC 2903 was estimated from the inclination corrected 
value by Jackson et al. (1991). The gas density for NGC 6946 is that listed by Young et 
al (1995) and is averaged over the NICMOS Camera 3 field of view (~ 45"). The gas densities 
for NGC 4449 are estimated using an average of values listed in Hunter & Thronson (1996). 
Gas densities for the LMC were taken from Luks & Rohlfs (1992) and Cohen et al. (1988). 
The molecular gas density for the NE position of M51 was taken from Kuno &: Nakai (1997). 
When necessary we have corrected published values so that we used the same conversion 
factor of N(H2) = 2.8 x lO^Ocm-^Ieo (K km s^^) (Scoville et al. 1987). 



